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Abstract 

The transverse energy spectrum in the unit rapidity window in pp collisions at ^/s = 
540 GeV is calculated to the next-to-leading order accuracy 0(a^) and compared to the ex- 
perimental data by UA(2) collaboration [[|. We show that the calculated spectrum starts 
matching experimental data only at relatively large transverse energy E± ~ 60 GeV and is 
in essential disagreement with it both in shape and magnitude at lower transverse energies. 
The data are well reproduced by HIJING Monte-Carlo generator indicating the crucial 
importance of all-order effects in perturbation theory as well as those of hadronization in 
describing the transverse energy production in hadron collisions at small and intermediate 
transverse energies. 
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In this note we calculate to the next-to-leading (NLO) order accuracy the transverse 
energy spectrum in the central rapidity window within perturbative QCD and compare it 
to experimental data obtained by UA2 collaboration 0. The NLO calculation of a generic 
jet cross section requires using a so-called jet defining algorithm specifying the resolution 
for the jet to be observed, for example, the angular size of the jet-defining cone, see e.g. 
|| , The cross section in question is calculated by integrating the differential one over the 
phase space, with the integration domain restricted by the jet characteristics fixed by 
the jet-defining alogorithm. Schematically the NLO distribution of the transverse energy 
produced into a given rapidity interval y a < y < y^ is given, to the 0(a^) order, by 

-^L = f D 2 PS 6(E ± - \P±i\°(ymin <Vi< Vmax)) 

dE± J d A pid 4 p 2 fr{ 

/do ^ 
D3ps ,4 , 4 , 4 S ( E ± ~ E \p±Mymin <y t < y max )) (i) 

where the first contribution corresponds to the two-particle final state and the second 
contribution to the three-particle one. The two-particle contribution has to be computed 
with one- loop corrections taken into account. 

In perturbative QCD one can rigorously compute only infrared safe quantities ||, in 
which the divergences originating from real and virtual gluon contributions cancel each 
other, so that adding very soft gluon does not change the answer. It is easy to convince 
oneself, that the transverse energy distribution into a given rapidity interval Eq. (1) 
satisfies the above requirement^ 

The calculation of transverse energy spectrum in pp collisions was performed in [|l| 
using the Monte-Carlo code developed by Kunzst and Soper 0, and a "jet" definition 
appropriate for transverse energy production Eq. (1). Recently the transverse energy 
production to the NLO accuracy was also discussed in ||, where the only the moments 
of transverse energy distribution were discussed and no comparison with experimentally 
observed transverse energy spectra was made. 

In Fig. [I] we compare the LO and LO+NLO transverse energy spectra for pp collisions, 
calculated as in |IJ, with the experimental data on transverse energy distribution in the 
central rapidity window \y\ < 1 and asimuthal coverage 7r/6 < ip < ll7r/6 at ^/s = 
540 GeV measured by UA2 Collaboration 0. We see that the perturbative LO+NLO 
calculations start merging with the experimental data only around quite a large scale E± ~ 
60 GeV. It is interesting to note, that it is precisely around this energy, that the space 

2 For a formal definition of infrare safety see, e.g., 
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Figure 1: Transverse energy spectrum in pp collisions calculated within LO+NLO accu- 
racy in perturbative QCD vs the experimental data by UA2 collaboration |2j 

of experimental events starts to be dominated by two-jet configurations 0. This means 
that only starting from these transverse energies the assumptions behind the perturbative 
calculation (collinear factorization at leading twist, explicit account for all contributions of 
a given order in a s ) are becoming adequate to the observed physical process of transverse 
energy production providing the required duality between the description of dominant 
configuration contributing to transverse energy production at this order in perturbation 
heory and the final state transverse energy carried by hadrons. At E± < 50 GeV the 
calculated spectrum is in radical disagreement with the experimental one both in shape 
and magnitude calculated and observed spectrum is very large indicating the inadequacy 
of the considered 0(a^) perturbative calculation in this domain. Let us mention here, 
that it is currently impossible to improve the results of the above calculation, because 
neither calculations of higher order nor infinite order ressummation for this process are 
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currently available []. 

In practical terms this means that additional model assumptions are needed to achieve 
agreement with experimental data strongly indicating that higher order corrections and 
higher twist effects have to be taken into account (in a necessarily model-dependent way) 
in order to describe them. In the popular Monte-Carlo generators such as PYTHIA ||] 
and HIJING jp such effects as multiple binary parton-parton collisions, initial and final 
state radiation and transverse energy production during hadronization are included. In 
Fig. ^|we compare the same experimental data by UA(2) with the spectrum calcuated 
with HIJING event generator. To show the relative importance of different dynamical 
mechanims, in Fig. |2| we plot the contributions from the hard parton scattering without 
initial and final state radiation, full partonic contribution and, finally, the transverse 
energy spectrum of final hadrons. We see, that taking into account additional partonic 
sources such as, e.g., initial and final state radiation, allows to reproduce the (exponential) 
form of the spectrum, but still not the magnitude. The remaining gap is filled in by 
soft contributions due to transverse energy production from decaying stretched hadronic 
strings. The results confirm those earlier communicated to the author by M. Gyulassy 
|IO|j . Let us note, that the spectrum calculated in HIJING is somewhat steeper than the 
experimental one. Additional fine-tuning can be achieved by probing different structure 
functions. 

The above results clearly demonstrate that in order to reproduce the experimentally 
observed transverse energy spectrum, one has to account for complicated mechanisms of 
parton production, such as initial and final state radiation accompanying hard parton- 
parton scattering, production of gluonic kinks by strings, as well as for nonperturbative 
transverse energy production at hadronization stage. This statement is a calorimetric 
analog of of the well-known importance of the minijet component in describing the tails 
of the multiplicity distributions, and |4| (see, however, the alternative explanation 
based on accounting for multipomeron contributions described in ||12|| ). 

Let us note that the result has straightforward implications for describing the early 
stages of heavy ion collisions. In most of the existing dynamical models of nucleus-nucleus 
collisions they are described as an incoherent superposition of nucleon-nucleon ones. As 
we have seen, to correctly describe the partonic configuration underlying the observed 
transverse energy flow in nucleon-nucleon collisions, mechanisms beyond conventional 
collinear factorization are necessary. In particular this indicates that the results obtained 
within minijet approach based on colinearly factorized QCD, see e.g. the recent review 
H| and references therein, must be taken with care. 

3 For discussion of the contribution of the initial and hnal state radiation see [B| 
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Figure 2: Transverse energy spectrum in pp collisions calculated in HIJING vs the exper- 
imental data by UA2 collaboration |2j 
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